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I. INTRODUCTION 



By taking Dirac's ideas [1] into account, R.P.Feynman explained how the non-relativistic 
Quantum Mechanics can be formulated from principles that make contact with the varia- 
tional principles of the Lagrangian Mechanics |§. Feynman showed that Quantum Mechanics 
can be based on the statement that the propagator, i.e. the probability amplitude of finding 
the system in the state q" at t", given that it was found in q' at t', can be obtained by 
means of the path integration: 

" i 



K(q" i"|q' t') = J Dq(t) exp 



h sHt)} 



(i) 



where S'[q(t)] is the functional action of the system. Since the path integral is a func- 
tional integration, one gives a meaning to Eq.(|I|) by replacing each path by a skele- 
tonized version where the path q(t) is represented by a set of interpolating points (q/c,£fc), 
k = 0, 1,...,N, q = q', q N = q", q fc = q(t k ). Then the functional action is replaced by 
a function S({(\k,tk}), and the functional integration reduces to integrate the variables q^, 
k — 1, N — 1 Q. Finally the limit Atk = tfc+i — tk — ► (i.e., N — > oo) is performed|. 
The function S({qk,tk}) is chosen to be P,f|: 



iV-l 



S({(lk,t k }) = J2 S (^k+itk+ihktk) (2) 



fc=0 



where S (q^k+itk+il^lktk) is the Hamilton principal function, i.e. the complete solution (in each 
argument) of the Hamilton- Jacobi equation that is equal to the functional action evaluated 
on the classical path joining its arguments. Thus the skeletonization (|2|) replaces each path 
q(t) by a succession of pieces defined by the system itself, which join the interpolating points. 
The skeletonized action (|^) retains the essential classical property of the functional action; 
namely it is stationary on the points interpolating the entire classical path between (q', t') 
and (q",t"). In fact , S({qk,tk}) is stationary when 



1 The convergence is assured by endowing the time with an imaginary part of proper sign. 

2 The rigorous mathematical meaning of this limit can be consulted in [fj] and references therein. 
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d d 
- — S(q k+1 t k+1 \q k t k ) + - — ^(qfctfelqfe-itfc-i) = 0, V/c (3) 

cqfc o<\k 
meaning that the q& are such that the final momentum of the classical piece between 
(qfc-i, tfc-i) and (q^, t k ), matches the initial momentum of the classical piece between (q^, t k ) 
and (qfc+i, ifc+i)- This continuity guarantees that the points {(q^,^)} are interpolating 
points of the entire classical path between (q',t') and (q",t"). 

Although a proper skeletonization for the path integral exists in the configuration space, 
the measure in Eq. ([!]) remains ambiguous. For instance, the finite propagator for a quadratic 
Lagrangian is known to be ||[| 



tf(q'V'|q'f 



det 



i 9 2 5(q' / r|q't / )M 1/2 / 



exp 



(4) 



2tx% dq"dq' 

This expression is also valid for the infinitesimal propagator of any classical system . The 
prefactor in Eq.([|) is the Van Vleck determinant ||, which takes part in the measure, and 
is nontrivial even for the short-time pieces of the skeletonization. 

A different kind of example is the (finite) Newton- Wigner propagator for the relativistic 
particle in flat space-time 

(in 1+1 dimensions), where S(q"t"\q't') = —mc(c 2 (t" — t') 2 — (q" — q') 2 ) 1 ^ 2 , and K\ is a 
modified Bessel function. In this case, the exponential of the Hamilton principal function 
does not cleanly appear in the propagator, and neither does it in the short time version 
(actually the propagator ([]) gets the form (|j ), not when t" — > t' but in the classical limit 
when the Compton wavelength h/(mc) goes to zero). Results of this sort could indicate a 
failure of (fj) to give the quantum propagator for an arbitrary system [fLOR . Anyway, it lays 
bare our complete ignorance of the measure in the representation ([!]). 

It was thought that a path integration in phase space could remedy this problem because 
there is a privileged measure in phase space: the Liouville measure c?q dp /{2nK) n (n is the 



dimension of the configuration space), which is invariant under canonical transformations. 
In this case, one should find an appropriate recipe for the skeletonization of the canonical 
functional action 

S[q(t), p(t)] = f (p(t) • q(t) - (q, p)) dt, (6) 
Jt' 

i.e., one should replace the functional in Eq.(^) by a function S^jq^, p&, tk}) of interpolating 
points for the path q(i), pit). The function ^({q/c, p^, tk}) should be dictated by the system 
itself. 



In Ref. [1 [| several recipes were essayed for a newtonian system moving on a riemannian 
manifold. Since the data (q^, p^; q^+i, Pk+i) overdeterminate the classical path between 
tk and tk+i, then the basic idea was to use the classical piece in the configuration space 
(just as in the previous case), together with the parallel transport of p(t&) = Pk along that 
classical piece. Of course, the parallel transport of p^ does not end in p^+i (unless the points 
interpolate the entire classical path between t' and t"). So the skeletonized path proves to 
be discontinuous in p (an unavoidable fact in phase space). The different recipes for the 
skeletonization came from the possibility of replacing the metric by a bitensor with the 
right coincidence limit. After the momenta were integrated on, an infinitesimal propagator 
similar to the one of Eq.(£|) was obtained. However the different skeletonizations reflected 
in a measure differing from the Van Vleck determinant by corrections associated with the 
curvature of the manifold. As a consequence, the Hamiltonian operator in the Schrodinger 
equation had a term proportional to h 2 R, where R is the curvature scalar (see also Ref. 



12|j for newtonian systems, and Ref. || for relativistic systems). 

The skeletonization proposed in Ref. [Tl[ successfully retains the covariance of the system, 



but it does not treat coordinates and momenta on an equal footing (a desirable feature in a 
canonical formalism). 



In Ref. [[L3|] the use of complete solutions of the Hamilton- Jacobi equation in the 



skeletonization has been suggested. A complete solution P^,[T5| 0(q, P,t) (the P's are 



n integration constants) can be regarded as the generator of a canonical transformation: 



p = <90/<9q, Q = d<p/dP, where (Q,P) is a set of classically conserved variables. Then 
dcf) = p ■ dq+ Q ■ dP — Hdt, and A0 P = 0(q", P, t") — 0(q', P, t') coincides with the canoni- 
cal functional action evaluated along a path such that P = const. Besides A0 P is stationary 
when P has the value corresponding to the classical path joining (q', t') and (q", £"); in that 
case A0 P turns out to be the Hamilton principal function JTSJ. These properties could make 
A0 P a candidate to take part in the skeletonization S({qk,Pk,tk})- 

But in order to obtain the propagator, some requirements concerned with the behavior 
at short times and the character of the substitution p — > P -which must be well defined in 
all phase space-, should be fulfilled by the complete solution to be chosen. In addition, the 
canonical coordinates and momenta should enter the skeletonization on an equal footing. 

The rest of the paper is devoted to emphasize the role played in phase space path 
integration by two related complete solutions of the Hamilton- Jacobi equation, which will be 
called Jacobi principal functions. In Section II a scheme of skeletonization putting canonical 
coordinates and momenta on an equal footing suggests the initial condition that must be 
fulfilled by the complete solutions to be used. In Section III the infinitesimal propagator 
induced by the path integration is obtained once the measure is worked up into a form that 
gives the same status to both arguments in the propagator. Section IV explains how to treat 
a classical system with arbitrary potential, in order to get the result (f§) for the infinitesimal 
propagator. Section V shows the ordering for the Hamiltonian operator that is induced by 
the propagator of Section III. The conclusions are displayed in Section VI. 

II. SKELETONIZATION IN PHASE SPACE: THE JACOBI PRINCIPAL 

FUNCTION 

In order to introduce a skeletonization procedure treating coordinates and momenta on 
an equal footing, one should define a recipe joining classical pieces determined by mixed 
boundaries (q^, Pk+i) or (p^, qjb+i). Then a path q(£), p(£) should be skeletonized by alter- 
nately giving the values of canonical coordinates and momenta at each tk, and replacing the 
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canonical functional action by something like 

(AT-2)/2 

S({<ik,Pk,tk}) = X! {J ( c l2k+2t2k+2\P2k+lt2k+l) + J (P2k+lhk+l\(l2kt2k)} ■ (7) 
k=0 

The building blocks J(q£'|p£) and J(p£'|q£) should be functions associated with the clas- 
sical system, making the skeletonized action stationary on the classical path in phase space. 
A comparison with Eq. (H) suggests that the stationary condition 
d d 

o J ( < i2k+2t2k+2\P2k+lt2k+l) + q </(p2fc+l*2fc+l \ ^khk) = 0, V/c (8) 

CP2fc+l CP2fc+l 

should mean that the final canonical coordinates of the classical piece between t 2 fc and £2*1+1 
coincides with the initial canonical coordinates of the classical piece between £2*1+1 and £2*1+2 • 
Once the stationary value for the momenta is replaced in Eq.(|7]), the skeletonization should 
go to the one of Eq.(^), so guaranteeing the continuity of both q and p at £2^+1 ■ 

If the system exhibits invariance under a general coordinate change, then the skele- 
tonization and the measure must preserve that invariance, in order that the quantization is 
independent of the chosen coordinates. Therefore both J's in Eq.(0) should be invariant. 



We are going to define dJ/dp to be the coordinate canonically conjugated to p; so 
we should look for coordinates transforming contravariantly to p. We will assume that 
the configuration space is a riemannian manifold Ai; thus, normal coordinates — which 
transform like the components of a vector at the origin of coordinates — could be introduced^. 
Let us choose a point O G M. as the origin of normal coordinates. Let {e a } be a basis for 
the tangent space Tq at O. To assign normal coordinates to a point P, consider the geodesic 



joining O and P and define a = s u, where s = J yj gijdq l dqi is the (invariant) length of 
the geodesic between O and P, and u eTq is the unitary vector tangent to the geodesic 



at O . The components a a (q J ) of the vector a G Tq are normal coordinates for P |16 



Actually, only a connection is needed to define normal coordinates. 
4 We assume a global topology such that any pair of points is joined by a unique geodesic. 
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By differentiating the invariant Hamilton principal function with respect to the normal 
coordinates of the (initial) final boundary, one gets the (initial) final momenta (—)p a = 
dS/da a ] thus p a are the components of a form p = p a e a e T Q . The set of canonically 
conjugated variables {(<J a , p a )} is invariant under changes qi — > q^'; (a a , p a ) only change 
under changes of the basis {e a } in T Q (and its dual basis {e a } in Tq). 

Now we will introduce two invariant Legendre transforms of the Hamilton principal 
function S(q" t"\q' f): 

W t"\ P > t>) = (S - ^ oA , (9) 

and 

Wlq'O = (s-^r«" a ) , (io) 

V uu / p'^=dS/da" a 

that will be called Jacobi principal functions. They can be regarded as the evaluation on the 
classical trajectory of a functional action that has been added with surface terms to make 
it stationary under variations with mixed boundaries left fixed. 

In spite of their definition in terms of the Legendre transform interchanging a a and p a , 
the Jacobi principal functions can be written, if preferred, as functions of a different set 
of canonical coordinates connected with the normal ones by means of a canonical point 
transformation 

= q 3 ( a a ) , Pj = — p a . 

The pj transform like the components of a form p G Tp on the coordinate basis . 

From the properties of the Legendre transform, one easily gets that J(q" t"\p' t') and 
— J(p" t"|q' t') generate contact transformations: 

„ d.J(^t"\p't') la dJ(q''t'\p't) 

p > = — ' = — w a — • ( } 

and 



a" a 



dJ(p"t"\q't') , dJ(p"t"\tjt' 



dp" a 3 dq'i 

When t" = t' they generate the identity: 

J(qt\pt)=p a a a (qi), (13) 

J(pt\qt) = -p a a a (qi) . (14) 



J(q" t"\p' t') and — J(p" t"|q' t') are complete solutions of the Hamilton- Jacobi equation 
in both arguments (take d/dt" and d/dt' in Eqs. (|9|) and (|TUD, and use Eqs. ( |TT]) and (|T2|) ) : 

^!iL0 = _ Jrtf ,W. O ME2 = E(,^L, ft ,fl (15) 

a t" a ^ at' o p a 

dJ(p t"\qt>) a dJ dJ(p t"\ q t>) dJ 

at" o p a of o q> 

By changing t" < > £' in Eqs. ([15]) and (0) one realizes that 

J(pt"|qt') = -J(qt'\pt"). (17) 

If the system is conservative, then the Jacobi principal functions depend on t" and t' 
only through the difference t" — t' . Therefore J(qt'\pt") = J(q, p, r = At), and the pieces 
of the skeletonized action (0) have the form 

J (q2fc+2t2fc+2|P2fc+lt2fc+l) +^(P2fc+lt2A:+l|q2fct2fc) 

= ^(q2fc+2, P2fc+1, T2fc+l) ~~ ^(q2fc, P2fc+1, — ^2k) (18) 



Thus the skeletonization (0) gets the form A0 P proposed in Ref. [|13|], although a specif! 
complete solution of the Hamilton- Jacobi equation is being used here. 



S 



Even for a non conservative system, the short time limit of the skeletonization (J7|) is 



^(cL2fc+2^2fc+2 |P2fe+1^2fc+l) + <^(P2fc+1^2A;+l| ( l2fc^2fc) 

— Pa2k+1 cr ' 1 2k+2 — H{ CT a 2 fc+2, Pj 2k+l ' ^2fc+l) At 2 fe+1 
~ Pa2fc+1 Cr£l 2fc + 2k, Pj 2k+l' ^2fc+l) At 2fc 



(19) 



On any smooth path it is valid that 0^+2 ~^ a 2k when i 2 fc+2 
skeletonized action goes to the canonical functional action: 



t 2 fc . Thus the 



■J ( c l2k+2't2k+2 |P2fc+1^2fe+l) + ^(P2fc+1^2fc+l | ^hk^k) 

> Pa2k+1 ( a 2k+2 - a 2k) -H( (T a 2fc ,p i2fc+1 ,t 2fc+ i) (t 2 k+2 ~ hk) ■ (20) 



The skeletonization scheme proposed in this Section is based on a pair of complete solu- 
tions of the Hamilton- Jacobi equation (in both arguments) that treat canonical coordinates 
a a and momenta p a on an equal footing; this fact is evident in the initial conditions ( |T5| ) 
and (PD- The Jacobi principal functions do not depend nor on the chart {qi} neither on 
the basis of the tangent space at O. They do depend on the way the configuration space has 
been cut from the phase space (the quantum propagation is not invariant under arbitrary 
canonical transformations). 



III. THE PROPAGATOR 



In order to give sense to the functional integration in phase space 
K{c[' t"\c[ t>) = J Vp(t) Vq(t) exp [1 S[q(t),p(t)] 



(21) 



our attention must now turn on the measure. Since the functional action is going to be 
replaced by an invariant skeletonized version, the "magical " measure T>p(t) T>q(t) should be 
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consequently replaced by a measure in the space of the variables {q2fc, P2fc+i}- This measure 
must be able to retain the geometrical behavior of the propagator, which is apparent in the 
manner of propagating the wavef unction: 

*(q", t") = J dq' K(qY'|qY) *(q', 0- (22) 

If the wave function \1/ is regarded as scalar, then the propagator should be invariant in its 
final argument but a density in its initial argument. However an scalar wavefunction would 
compel us to use an invariant measure /x(q) <iq in the inner product of the Hilbert space 
(the density fi would be ultimately dictated by the result of the path integration ||). So it 
may be more convenient to regard the wave function as a density of weight 1/2. In this case 
the inner product of the Hilbert space is 

= J dc^m* $, (23) 

whatever the generalized coordinates describing the system are. Thus the propagator in 
Eq.(p2"D must be a density of weight 1/2 in both arguments. 

The issue of the measure can be studied at the level of an infinitesimal propagator. In 
fact, due to the composition law 

W t"U = / #(q"t"| qi v-i*iv-i) rfqjv-i if(qAr-itAT-i|qAr-2*7V-2) 

dq 2 tf(q 2 * 2 |qiti) dqi #(qi*i|qV) (24) 

— which holds whenever the added paths go forward in time — , the finite propagator can 
be retrieved by composing infinitesimal propagators. If t" — t' = e is infinitesimal, then 
one should only integrate p at some intermediate time t. However the measure dp" does 
not allow for a propagator behaving like a density in its arguments q' and q" . The use of 
d n p%' instead of dp^ is suitable when the wave function is regarded as a 
scalar, because the propagator will result a density in q' . However, if the propagator has 
to be a density of weight 1/2 in both arguments, then the Jacobian in the previous measure 
must be splitted in two factors that will give an equal weight to p' and p" : 
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cZV- 



det 



dp', 



dp'a 



(*) 



det 



dp'- 



1/2 



d pa 



(t) 



d n P a t] 



det 



dp'j 



1/2 



det 



'd 2 J(q"t"\pt)' 
dq"i dp a 



1/2 



d n Pa t] 



det 



d 2 J(pt\g't') 
dq' i dp a 



1/2 



(25) 



Concretely, the infinitesimal propagator has the form: 



K(q" *" = f + e | q' = J 



d n Pa 



det 



a 2 J(q"t"|P*) 

dq" i dp a 



1/2 



det 



d 2 J(pt\q't') 

dp a dq' i 



1/2 



exp 



- (J(q"t"\pt) + J(pt\4t')) 



(26) 



where t is prescribed to be the mid time: t = t' + (e/2) — t" — (e /2). When e = 0, one gets 
the orthonormality relation between eigenstates of the operator q (see Eqs.(O) and (p^)). 



Since canonical coordinates and momenta were treated on an equal footing, one realizes 
that the propagator in the p a - representation which is consistent with Eq.(^) is 



K{p" t" = t' + e\ p' t') = J 



d n q j 



det 



d 2 J(p"t"\qt) 
dp 'L 9q j 



1/2 



det 



d 2 J(gt\p't') 

Qqd Qp' a 



1/2 



cxp 



(J(p"t"\qt) + J(qt\p't')) 



(27) 



It will become clear in Section V that the prescription of mid time in Eqs.( p6|),(p7|) 
together with the splitting of the Jacobian, guarantee the hermiticity of the Hamiltonian 
operator and the unitarity of the evolution. 

IV. CLASSICAL SYSTEMS 

To illustrate the use of Eq.(|6|), let us consider a one-dimensional classical system gov- 
erned by the Hamiltonian 
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B = £ + ffe). (28) 



We will show how to manage the integration in Eq. (|2q) in order to get the infinitesimal 
propagator in the form of Eq. (|4|) . 

Since the metric in the Hamiltonian is a standard euclidean metric (g^ = 5^), the 
coordinate q is the normal coordinate. The Jacobi principal function J(q, p, r) can be guessed 
by writing 

v 2 00 

J(q,p,r) = pq - £-t + ^2Ji(q,p)r l . (29) 
Then one solves the Hamilton- Jacobi equation order by order in r, and obtains 

Jl = ~V(q) 
p V'(q) 



2m 

p 2 V"(q) V'(q) 2 



6m 2 6m 



5 ' 120 m 4 W 15 m 3 40 m 3 15 m 2 ' 1 j 



The recurrence formulae is 

2(Z + l)m ^ ag 

where Jo = pq. Each J is polynomical in p. Let us concentrate on the higher degree 
contributions; their addition is 

We are going to replace this result, and d 2 J/dqdp = 1 + 0(e 2 ), in the integrand of 
Eq. (p6|) . After the substitution p — > P = pe/(2m), one gets 
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K(q" t" = t' + e\q' t') 
dP exp 



m 
irhe 

2 



%2m / Aq\ 



exp 



16 

2hP 



im(Ag)" 
2fte 



/'/ rq'+P \ 

/ V(q)dq+ V(q)dq) + 

Ja"-P Jo' 



+ 0{e 2 ) . (33) 



The contribution of the potential to the phase will be expanded about P = Aq/2: 
1 / ri" ri'+ p \ 

= V — AV (P-**) + -^-AV (p- + 

Aq \ 2 j + (Ag) 2 V 2 ) ^ ' 



(34) 



where 



* 3 i // V(5) d5 ' 



AV = V - V 



g' + q" 



(35) 



Those contributions that were not explicitly written in Eq.fl3"3|) can be controlled by 



means of the result 17 



J dx x 21 exp 



i2m 



he 



-x 



\2mj 



7+ \ 



(36) 



Then the leading contribution to the integration (P6|) is 

i / m(Aq) 2 



K(q" t" = t' + e\q' t' 



m 



2i7ihe 



exp 



h \ 2e 



(37) 



The infinitesimal propagator ( [371) nas the form (f|). In fact the phase in Eq. (|37|) solves 
the Hamilton- Jacobi equation in each argument at order e (for all values of q' and q")f\- The 
Schrodinger equation is satisfied at the lowest order in he/m. 



5 Although the phase in Eq.(|37|) has the merit of being a complete solution of the Hamilton- Jacobi 
equation at order e -it is the Hamilton principal function at that order-, the integration on q in 
the composition of infinitesimal propagators (Eq.(f24j)) will be not sensitive to a replacement of V 
by V((q" + q')/2), or (V(q') + V(q"))/2, etc. as a consequence of the result @). 
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V. OPERATOR ORDERING 



Each recipe to path integrate implies an operator ordering for the Hamiltonian in the 
wave equation. Our interest in this Section is to find the operator ordering associated with 
the infinitesimal propagator (|26| ). Let us derive Eq.(p6|) with respect to e, at e = : 



i h ^-K(q" t' + e\q r t' 
oe 



dp 



exp 



p Aa 



ik ( d 2 H(q",p) d 2 H(q',p) 



+ 



[4 \ dq"dp dq'dp 
where the short-time approximation (|19|) has been used 



+ -(H(q",p) + H(q',p)) 



(38) 



Eq . (|38|) is linear in H. If H can be expanded in a power series, then it will sufficient to 
handle the ordering for a Hamiltonian 



H = q m p k . 



(39) 



If q is a normal coordinate, then 



i h ^-K(q" t' + e\q't' 
Oe 



e=0 



dp 
2nh 



exp 



p Aq 



ih k m k _ x 



p k-l (_ q " m-l + g > m-l) + |_( ? "m + g /m) 



(40) 



Taking into account the Eq.(P2| 



n p Aq 



ih k m ^ k _ 1 f j, m _ 1 | j m _x\ L P_ ^// m _|_ ^ m -A \gfq' 



pfc-l (_ g " m-l + q , m-l 



n m-iih k 171 ( h d \ k 1 [ dq' dp 



-q 



4 1 i d l" 



2nh 



exp 



p Aq 
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+ 



ih k m ( h d 



fc-i 



i dq" 



dq' dp , m _! 
-TT- 9 exp 

ZTlh 



- p Aq 



2 q 



h d \ f dq' dp 



i dq" J J 2nh 



exp 



-p Aq 



1 fh d \ f dq' dp , m 

+ 2Uw)j^ q exp 



-p Aq 



ih k m 



+ ^(P h q m + q m p k ) W,0). (41) 



This means that the Hamiltonian operator is 

1 / .™ „™ ih k m 



H = ~{p k q m + q m p k ) + | 



(42) 



The operator H is Hermitian thanks to the mid time prescription in Section III, which 



gave an equal weight to the terms depending on q' and q" in Eq . (|38|) . 



VI. CONCLUSIONS 

We have proposed a scheme to path integrate in phase space, which is applicable to 
Hamiltonian systems whose configuration space is a manifold where normal coordinates 
(i.e., coordinates behaving like the components of a vector in the tangent space at the 
origin) can be introduced. The skeletonization is based on the invariant Jacobi principal 
functions — those related with the variational principles of Mechanics for mixed boundaries 
left fixed — , and the measure gives to the propagator the character of a density of weight 
1/2 in each argument. The so obtained infinitesimal propagators (^) and ( p7| ) naturally 
satisfy the Schrodinger equation, once the Hamiltonian operator is build in agreement with 
the operator ordering (|42"D induced by the path integral recipe. 

The infinitesimal propagators ( p6| ) and (^) can be read in terms of the modes 

F P (q j ,t) = 

15 



det 



1 dV^p,^ 



2nh dqi dp a 



exp 



-J(q,p,t) 



(43) 



which are well behaved on all phase space at short times (the matrix d 2 J j dq 1 dp a is not 
singular because the relation between initial and final momenta is biunivoque at short times). 
At t = these modes are eigenfunctions of the momenta p a = —ihd/da a based at the origin 
O, because of the boundary condition J(q, p, t = 0) = p a c a . Therefore J r p (q : ',t) is a short- 
time approximation for < g J \U(t)\p a >, and {F p } is a basis of short-time solutions of the 
Schrodinger equation whatever the Hamiltonian system is. A change of the origin O implies 
a change of the basis {.Fp}; of course, all bases {Fp} are equally good for expanding the 
propagator. 

For free systems it is J(q, p,t) = p a a a — H{p a )t , and the modes ([43]) are certainly 
exact solutions of the Schrodinger equation. They span the basis of eigenstates of the 
(conserved) momenta p a . In this case the propagator (^) is exact (i.e., it is the finite 
propagator). In particular, the Newton- Wigner propagator (|5|) for the relativistic particle 
-H(p) = (p 2 + m 2 ) 1 / 2 - can be obtained by integrating on the momenta in Eq.(^) ||. 

In the case of the classical system of Section IV, the modes T v satisfy the equation (use 
Eq.©) 



' d_ tf_&__ v{ A T ( o 2 jy 1/2 

dt 2m dq 2 I p 2m \ dqdp J dq 2 



(—) 

\dqdp J 



1/2- 



T p . (44) 



This equation is typical for any phase J being a solution of the Hamilton- Jacobi equation, 
and is commonly used to remark the semiclassical character (h — > 0) of wavefunctions having 
the form (^3|) . However, as was already stated, the Jacobi principal function J(q, p, t) confers 
an additional property to the modes (f43|) -which is the one exploited in this paper- and allows 
a different reading of Eq. (|44]) : the modes T p are short-time solutions of the Schrodinger 
equation, for any value offi. In fact, the results in Section IV show us that the rhs in Eq.(f44D 
is zero for a quadratic potential^ [T v is an exact solution), h 2 t 4 /(Sm 3 )(V"') 2 J r p + 0(t 5 ) for 



e HV(q) = mi0 2 q 2 /2, the Jacobi principal function is J(q,p,t) = —p 2 t&n(u;t)/(2mu})+qp/cosLut — 
mujq 2 tan(cjt)/2. Thus J(q,p,t) reflects the equivalent roles of q and p in both the initial condition 
and the Hamiltonian. 
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a cubic potential, and h 2 t 2 /(8 m 2 )V""J r p + 0(t 3 ) in a more general case. 



The substitution of the functional action by a skeletonized version in the discrete-time 
approximation is one of the ways to give a meaning to the functional integration (ED). A 



different approach to the same problem is the operator symbol method [0 , where the path 
integral results from the product of the symbols associated with the short-time evolution 
operator. The product of symbols involves an integration containing the information about 
the operator ordering, which amounts to the prescription of the skeletonization in phase 
space. The different rules to generate the ordering for the quantum operator g associated 
with a function g(q,p) in phase space, can be summarized as follows fTT 



9 



f (-i— -i— 1 exp (-—-?—) 
\ dq ' dp J \ 2 dqdp ) 



g(q,p), 



(45) 



where g in Eq. ([45]) is the normal form of the operator (the power series expansion where 
the q precede the p), and f(u,v) contains the information about the ordering. Some well 
known rules of ordering are ||20|| : 



Name 




g(q,p) = q m p k 


Standard 


exp[^huv] 




Anti-Standard 


exp[— Tjhuv] 




Symmetric 


cos[^huv] 




Weyl 


1 




Born- Jordan 


2 (huv)' 1 sin[^huv] 


k 

_1_ £ pk-l g mpl 



The Weyl ordering is equivalent to a skeletonization where the Hamiltonian is evaluated 
in (q" + q')/2. The symmetric ordering corresponds to replace the Hamiltonian by [H(q",p) + 
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H(q',p)]/2 PH| . In this sense, the skeletonization prescribed in this paper seems to be related 
with the symmetrization rule (see Eg. (|T9|) ) . Efectively, the ordering (^) for the normal 
coordinates and their conjugated momenta begins with a symmetrized contribution coming 
from the mean Hamiltonian in Eq.(|3^). However, our prescription includes a non-trivial 
measure (the Jacobians in Eq.(|25D), which is needed in order that the wavefunction retains 
its condition of being a density of weight 1/2. This measure is responsible for the second 
term in the ordering (fl2f). Thus the function f(u,v) associated with the ordering (|4~2"D is 

111 

f(u,v) = cos[-huv] + -huv sm[-huv]. (46) 

jU Zi £ 

Naturally / fulfills the requirements 

lim / = 1 lim / = (47) 
h^o h—>o 



that guarantee the classical limit [19 . 
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